Director precession and nonlinear waves in nematic liquid crystals under elliptic shear 
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Elliptic shear applied to a homeotropically oriented nematic above the electric bend Freedericks 
transition (FT) generates slow precession of the director. The character of the accompanying non- 
linear waves changes from diffusive phase waves to dispersive ones exhibiting spirals and spatio- 
temporal chaos as the FT is approached from above. An exact solution of the flow alignment 
equations captures the observed precession and predicts its reversal for non-flow aligning materials. 
The FT transforms into a Hopf bifurcation opening the way to understand the wave phenomena. 



Liquid crystals (LCs) exhibit a multitude of interest- 
ing nonlinear dynamical phenomena, like electrically and 
thermally driven convection, flow-induced and optical in- 
stabilities 0J^. Nonlinear dispersive waves as observed 
in particular in excitable and oscillatory chemical reac- 
tions ^ are rarely seen in LCs. An exception are the 
2D waves that modulate the oscillating bimodal pattern 
found in electroconvection of nematics Q for which only 
a phenomenological description has been developed. 

Actually, nonlinear dispersive waves are obtained in 
a much simpler situation when a local precession of 
the director is generated by an elliptic shear flow in a 
homeotropically anchored, Freedericks distorted nematic 
slab. As will be shown here, the dynamics can then be 
nearly conservative, in contrast to the well studied situ- 
ation where the director is set into motion by a rotating 
magnetic or electric field and the dynamic phenom- 
ena are limited by the strong dissipation. In the presence 
of elliptic shear the FT transforms into a Hopf bifurca- 
tion and the observed transition to chaos can possibly be 
related to the Benjamin- Feir instability, see e.g. 

Elliptic shear flow in homeotropically aligned nematics 
has been studied intensively in the past in view of an in- 
stability leading to roll patterns For that purpose 
rather large shear angles are needed and no external field 
is applied. The only previous study of the Freedericks dis- 
torted state in the presence of elliptic shear was, to our 
knowledge, carried out by Dreyfus and Pieranski. They 
observed some of the phenomena discussed below with- 
out a conclusive interpretation |^ . 

In this work the nematic layer (thickness d) has been 
sandwiched between clean, Sn02 coated float glass plates 
(xy plane) resulting in homeotropic anchoring (nematic 
director fi perpendicular to the bounding plates) of the 
substance Phase 5A (Merck). A voltage U applied 
across the sample induces at a critical value Up a bend 
Freedericks transition (FT) (the dielectric anisotropy Ea 
of Phase 5A is negative ) . Above U f the director tilts 



away from the z direction. The tilt angle 9 can be con- 
trolled by varying U. The degeneracy with respect to 
the azimuthal angle leads to the well known "schlieren 



texture" seen between crossed polars |10| . It contains de- 
fects {umbilics) with topological charge ±1 depending on 
whether the in-plane director makes a 27r or — 27r rotation 
on a closed loop around them. These can be distinguished 
by rotating the crossed polars. 

Elliptic shear flow has been generated by oscillating the 
bottom glass plate in the x direction [x{t) = Ax cos(wt)] 
and the top plate along y [y{t) — Ay cos{ujt + $); unless 
stated otherwise the phase shift will be <& = — 7r/2]. The 
oscillations were produced by loud speakers. The preces- 
sion of the director was observed as the local oscillation 
of the transmitted light intensity that results in wave-like 
propagation of bright and dark domains. The recordings 
were studied by digital image analysis. Details of the 
experimental setup are described in 

At sufficiently high voltages the director orientation 
varies slowly in space and precesses almost homoge- 
neously in time (Fig.|l| bottom). Around 1.2Uf inho- 
mogeneities emit traveling waves and umbilics generate 
spiral waves, very similar to those observed in oscillatory 
and excitable chemical reactions ||^. The longer waves 
in the background originate from the lateral cell bound- 
aries (Figj^ middle). At lower voltages spiral pairs seem 
to be created spontaneously (without umbilics) and one 
observes spatio-temporal chaos (Fig.|l| top). 

In Figjl the dependence of the precession frequency fl 
on the voltage is shown for two temperatures. Clearly 
17 behaves differently in the different regimes (see FigJ^) . 
The dependence of ft on the frequency and amplitude of 
the applied shear is given in Fig. ^. The linear dependen- 
cies can be understood from symmetry arguments and 
dimensional analysis as shown in the following. 

In order to describe the director precession we look for 
solutions of the nematodynamic equations where the di- 
rector n and the velocity v — {vx,Vy,0) depend only on 
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the coordinate z and time t. Measuring lengths in units 
of d, times in units of we introduce 

ki Ki^jKzz , A = azlo-2 , = 71^^/7^33 , 71 = 03 - 02 
e = l/(rdt^) , = {72/;7| , C/f = V^33^7(eo|ea|) 

with iCii the elastic constants and a; the Leslie coeffi- 
cients. Setting N = ^ iuy , V = Vx + iVy, one can 
write the nematodynamic equations concisely in complex 
notation. The torque balance equation is 

nz ( 1 + A,— — A 

-e {kiUzUz^zzN - [Nzz ~ N{NN,,,+NN^,,)/2] 
+ (1 - k2)[N{B^ + iB,z)/2 + iN.zB] - TT^e^nlN} , (2) 

with Uz = y/1 - |iV|2, B = i{NN,z - iViV,z), and de- 
note derivatives. 

For frequencies such that the viscous penetration depth 
satisfies I — ^/yi/puj » d [f — uj/{2tt) below about 
40 kHz for d — 2Q /im], and small dimensionless shear 
amplitudes a = A^jd one has a linear flow field, i.e. 
V-iz — a(cost -|- i5sint) where h — Ay/A^ is a measure 
of ellipticity {b = 1 for circular shear) . 

It is useful to first neglect any space dependence, which 
is a good approximation when the Ericksen number a/e 
^ is large. Consistently, one then also has to discard 
the electric field, so that in Eq.(^ one is left with the 
terms in the first line. Introducing angles by writing 
N = sin 9 exp i(f) one can rewrite the equation as 

6,t = a' (cos^ 6* — A sin'^ 0) [cost cos -f 6 sin fszn0] , 
(j),t = a' cot 0[— cost sint/) -f bsint cos (/)] , (3) 

where a' — a/(l — A). For rectilinear shear (6 — 0) one 
recovers the flow-aligned solution cot^ 9 — X, cj) — 0. 

Equations (|^) represent a conservative, reversible dy- 
namical system. For A = the director is advected pas- 
sively by the velocity fleld and Eqs. (j^) separate into 

dt (tan 9 cos (f>) — a cos t , dt (tan 6 sin (p) ~ ab sin t . (4) 

The solutions for arbitrary initial conditions can easily 
be written. They describe simple, closed 27r— periodic 
orbits, which either include or exclude the origin 
(homeotropic orientation). Clearly this case can be gen- 
eralised to arbitrary time dependence of the flow. 

For circular flow (6=1) Eqs.(^ are integrable even 
for A 7^ 0. Then the terms in square brackets become 
cos (0 — t) and — sin {(f> — t), respectively. Introducing 
the phase lag (p ^ <f> — t the equations become au- 
tonomous. Transforming them into the 2nd-order ODE 

Lp^tt — {'~P,t +l)(2<y£',t -f 1) cot ip — Xa"^ simp cos Lp (5) 

one can verify that the quantity 



(^.«+l-A«-7in^^)^ (6) 
(2<p,t — Aa'^ sin p) sin p 

is a constant of motion. Solving for Lp,t one obtains the 
2^ period T of the motion as an integral which can be solved 
analytically giving 

T= dip / if, t = 2t: / 1 - Xa'^ (7) 
Jo 

(independent of C!). For A 7^ the orbits are in general 
quasi-periodic. Thus, in addition to the rapid oscillation, 
the director performs a slow precession with frequency 

— (1 — 27r/T)a; « {X/2)a'^LU in physical units. The pre- 
cession is for flow-aligning materials positive (same sense 
of rotation as the elliptic shear) and negative otherwise. 
In Fig. ^ some typical orbits are shown. 

Simulations show that for b ^ 1 (not too small) and 
A 7^ the slow precession persists and has the same sign 
as for b = 1. One has ~ {X/2)a''^bu) , which can also be 
derived analytically for small shear flow amplitudes. 

The dependence of 17 on the ellipticity is explored 
experimentally most easily by keeping A^ — Ay fixed 
and changing the phase shift $. This driving is actu- 
ally equivalent to one with phase shift — 7r/2 and A'^ = 
cos($/2), A'y = \/2Aj, sin(<I>/2), as can be seen 
by transforming into a rotated coordinate system. Thus 
a^b = {AxAy/d'^)\sin<^\. In Fig. || (insert) the experi- 
mental results are compared with the theoretical predic- 
tion, ri/cij was always found to be positive as expected 
for the flow-aligning case. From the magnitude of at 
U/Uf « 2 and T = 25 . . . 30°C in Figs. 2, 3a, and 3b. 
we find A — 0.008 ±0.004, which appears consistent with 
the very small values found in and falls within the 
experimental uncertainty of those measurements. 

Now we discuss the effect of the terms in curly brack- 
ets in Eq. (|^). These terms represent singular perturba- 
tions. They introduce dissipation into the system, which, 
together with the boundary conditions, produces the at- 
tractors. For 6 = 1, a << 1 and at lowest order (in e), 
this reduces to a (z-dependent) selection from the family 
of solutions of the conservative system. Although this 
approximation breaks down near the bounding plates 
it can give useful results (see the rectilinear case [Q). 
One then expects that (aside from the slow precession) 
the director performs small-amplitude oscillations around 
the director distribution 9q{z) obtained from equating 
the curly brackets in Eq.(^ to zero (with homeotropic 
boundary conditions iV = 0), which corresponds to the 
usual (static) bend Freedericks distortion. In this ap- 
proximation the frequency of the slow precession should 
be as calculated above. We also analysed slow modu- 
lations in the xy plane and found them to give rise to 
diffusive phase waves as observed well above Up ■ 

The experimental results indicate that corrections to 
the above behavior are important, in particular near 
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the FT. We have analysed corrections only in limiting 
cases and found the contribution to the slow precession 
to be positive. In particular, one should note that in 
the presence of circular shear the FT is transformed into 
a Hopf bifurcation, as can be seen purely from symme- 
try arguments: below the transition the director rotates 
symmetrically around the homeotropic orientation. At 
the FT this symmetry is broken and generically a new 
frequency develops (no chiral symmetry in the rotating 
frame) which leads to the precession. 

Consequently, near the FT a description in terms of a 
complex Ginzburg-Landau equation should be possible. 
For low frequencies (e ~ 1) and near the FT we have 
performed a fairly full analysis. We find the FT to be 
slightly delayed and modulations in the xy plane to ex- 
hibit strong dispersion, which presumably leads to the 
observed modulational (Benjamin-Feir) instability. For 
l>6>5c~0.5 the behavior remains qualitatively simi- 
lar. At be there is a crossover from oscillatory to excitable 
behavior of the system . 

In summary: the theory describes the slow precession, 
at large fields essentially quantitatively. The scenario 
of the waves on the background of the slow precession 
(diffusive phase waves at large fields changing to am- 
plitude waves with dispersion and eventually Benjamin- 
Feir chaos at low fields) can be understood qualitatively. 
There remains to be done a quantitative analysis at low 
fields as well as an experimental test of the most provoca- 
tive prediction, namely the reversal of the slow precession 
for non-flow aligning materials in situations where the 
elasticity-induced effects are small. 

We have observed similar behavior in the nematics 
MBBA and DOBHOP (e^ < 0). Also, in 5CB (£„ > 0) 
we generated the director precession by inducing the di- 
rector tilt by surface effects: the bottom plate of the hor- 
izontal cell was heated above the phase transition point 
(34.5° C) so that a nematic-isotropic interface developed 
in the xy plane. The surface energy of this interface is 
minimized when the director encloses an angle ~ 64° with 
the surface normal 

We have also performed experiments with linear me- 
chanical vibration along z (compression) generated by 
piezo crystals (/ « 5 - 100 kHz) attached to one of 
the bounding plates. This presumably induces Poiseuille 
flow. As before, the slow precession occurred only in the 
Freedericks distorted state. Here the phase waves are 
typically emitted from certain locations in the form of 
target patterns, which presumably result from spatial in- 
homogeneities in the flow. The waves behave diffusively 
(even near Up), which is probably due to the fact that at 
the high frequencies used the elastic contributions to the 
precession are irrelevant For more details, see jl^ . 

Director precession and phase waves have been ob- 
served previously in cells that were excited piezo electri- 
cally at frequencies around 50 kHz . The piezo crystal 
formed one of the bounding plates. Phenomena reminis- 



cent of the phase waves were also seen in planar and 
homeotropic cells without electric field at frequencies 10 
kHz < / < 1 MHz ||l^ . There the waves originated from 
orientational defects at the surface. Although in these 
cases the precise form of the excitation was not clear we 
suggest that the mechanism presented here forms the ba- 
sis of the phenomena. 

Finally we note that the precession should be observ- 
able also in smectic C layers excited elliptically. Possibly 
the effect can be used to measure the analog of a3/a2- 
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FIG. 1. Snapshots in the different regimes (crossed po- 
lars). Circular shear. Ax = Ay = 3.4 fim, f = 155 Hz, and 
T = 28.2°C. 
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FIG. 2. The precession frequency f2 as a function of the 
applied vohage for T = 29°C, Uf = 8.4 V (•) and T = 53°C, 
Uf = 10.7 V (o). The other parameters were: d = 20 ^m, 
Ax = 3.1 fim, Ay = 3.5 fim, f = 155 Hz and $ = 80°. In 
the insert we show as a function of the phase shift $ for 
Ax = Ay = 4 lira, f = 122 Hz, d = 20 ^lm, U/Uf = 2.3 and 
T = 24.5°C (*). The continuous line shows VL = 0.102| sin<l>|. 
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FIG. 3. Precession frequency vs. driving frequency and 
amplitude. The parameters are: a.) U/Uf = 2.3, d = 20 
/im, Ax = 0.8 yLm, Ay = 1 /im and T = 25.5°C. b.) Circular 
shear, d = 20 /^m, / = 122 Hz, T = 28.5°C and U/Uf = 2.3. 




FIG. 4. Three orbits of the in-plane director under circu- 
lar flow for A = 0.2, a — 0.5. The thick circle represents a 
pure rapid rotation around the z axis relevant below the FT. 
(In the elliptic case this orbit transforms into an approximate 
ellipse.) Other orbits exhibit the slow precession. We show 
examples with a small average tilt (dotted), expected to be 
relevant slightly above the FT, and with a large tilt (dashed). 
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